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Abstract 

The bicycle model of Lorentz noninvariant neutrino oscillations without neutrino masses nat- 
urally predicts maximal mixing and a 1/E dependence of the oscillation argument for v ^ — > v T 
oscillations of atmospheric and long-baseline neutrinos, but cannot also simultaneously fit the 
data for solar neutrinos and KamLAND. Within the Standard Model Extension, we examine all 
19 possible structures of the effective Hamiltonian for Lorentz noninvariant oscillations of mass- 
less neutrinos that naturally have a l/£ dependence at high neutrino energy. Due to the lack 
of any evidence for direction dependence, we consider only direction-independent oscillations. 
Although we find a number of models with a l/£ dependence for atmospheric and long-baseline 
neutrinos, none can also simultaneously fit solar and KamLAND data. 
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1 Introduction 



Neutrino data from atmospheric, long-baseline, solar and reactor experiments are easily explained 
by oscillations of three active, massive neutrinos pQ. Lorentz-invariance and CPT violating inter- 
actions originating at the Planck scale can also lead to neutrino oscillations. The Standard Model 
Extension (SME) [2J includes all such interactions that may arise from spontaneous symmetry 
breaking but still preserve Standard Model gauge invariance and power-counting renormalizability. 
Studies of neutrino oscillations with Lorentz invariance violation have been made both for mas- 
sive IH [S] and massless [SI El [5] neutrinos. A model with nonrenormalizable Lorentz invariance 
violating interactions and neutrino mass has also been proposed [9J. However, no viable model 
has been found that does not require at least one nonzero neutrino mass. The purpose of this 
paper is to determine if Lorentz invariance violation alone can account for the verified oscillation 
phenomena seen in atmospheric, long-baseline, solar and reactor neutrinos. We do not attempt to 
fit the possible oscillation signals seen in the LSND [TU] and MiniBooNE [11J experiments. 

In the SME, the evolution of massless neutrinos in vacuum may be described by the effective 
Hamiltonian [6] 

(heff)ij = E5ij + - [al Ptl - d£p,& v ] .. , (1) 

where p^ = (E,—Ep) is the neutrino four-momentum, p is the neutrino direction, i, j are flavor 
indices, and — > —cll for antineutrinos. The coefficients ai have dimensions of energy and the cl 
are dimensionless. Direction dependence of the neutrino evolution enters via the space components 
of ai and cl, \i or v = X, Y, Z, while direction independent terms have p, = v = T '. The Kronecker 
delta term on the right-hand side of Eq. JIJ) may be ignored since oscillations are insensitive to 
terms in h e ff proportional to the identity. 

The two-parameter bicycle model [6] can be defined as follows: {ci)ij has only one nonzero 
element in flavor space and the only nonzero {a£)ij are (ai) e/1 = (az,) er . These interactions can 
be nonisotropic, which could lead to different oscillation parameters for neutrinos propagating in 
different directions. In Ref. [8] it was shown that the pure direction-dependent bicycle model is 
ruled out by solar neutrino data alone, while a combination of atmospheric, solar and long-baseline 
neutrino data excludes the pure direction-independent case. A mixture of direction-dependent 
and direction-independent terms (with 5 parameters) is also excluded when KamLAND data are 
added 0. 

The key feature of the bicycle model is that even though the terms in h e ff are either constant 
or proportional to neutrino energy, at high neutrino energies there is a seesaw type mechanism that 
leads to 1/E behavior for the oscillation argument for atmospheric and long-baseline neutrinos. In 
this paper we examine the general case of direction-independent Lorentz invariance violation in the 
Standard Model Extension for three neutrinos without neutrino mass, i.e., Eq. (pQ) with only c^ T 
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and aj^ terms. We do not consider possible direction-dependent terms since there is no evidence 
for direction dependence in neutrino oscillation experiments (see, e.g., the experiments in Ref. |12j 
and the analysis of Ref. [B]). For notational simplicity we henceforth drop the L subscript and T 
superscripts from the cj 7 * and a 7 in our formulae. 

We first look for textures of the Cj,- in flavor space that allow &1/E dependence of the oscillation 
argument at high neutrino energy. We then check the phenomenology for atmopheric, long-baseline, 
solar and reactor neutrino experiments. We were unable to find any texture of h e ff that could 
simultaneously fit all the data. 

In Sec. 2 we review the constraints on the direction-independent bicycle model. In Sec. 3 
we list all possible textures of the c coefficients and find which ones allow a 1/E dependence of 
the oscillation argument at high neutrino energies. For those that do, we first check the oscillation 
amplitude for atmospheric and long-baseline neutrinos, and if suitable parameters are found we then 
check the ability of the model to fit KamLAND and solar neutrino data. In Sec. 4 we summarize 
our results. 



2 Neutrino oscillations in the bicycle model 



As an illustrative analysis, we begin with a review of the direction-independent bicycle model and 
show how it is inconsistent with a combination of atmospheric, long-baseline and solar neutrino 
data. 

Neutrino oscillations occur due to eigenenergy differences in h e f f and the fact that the neutrino 
flavor eigenstates are not eigenstates of h e ff. In our generalization of the direction-independent 
bicycle model, 

/-2cE + 2a u a 12 a 13 \ 



eff 



V 



(2) 



ai2 

ai 3 / 

where the c term is CPT-even and the a,ij terms are CPT-odd. The simple two-parameter bicycle 
model [6] has ai3 = a\ 2 and an = 0. We allow a\2 to be different from 013 so that mixing of 
atmospheric neutrinos may be (slightly) nonmaximal. The an term allows an adjustment of the 
oscillation probabilities of low-energy solar neutrinos [6] . 

For this h e ff there are two independent eigenenergy differences Ajk = Ej — Ej- given by 



A 2 i = V(an - cE) 2 + a 2 + cE- a n , A 32 = ^(a u - cE) 2 + a 2 - cE + a n , (3) 
where a = \J a'f 2 + a'f 3 . The effective Hamiltonian is diagonalized via U T h eff U by the energy- 
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dependent mixing matrix 



U 



( — cos 9 
sin <j) sin ( 
cos 6 sin ( 





COS (j) 

— sin< 




(4) 



where 



sin 



tan i 



1 
2 
Ql2 
«13 



1 + 



an — cP 



^(an - cE) 2 + a 2 



The off-diagonal oscillation probabilities are 



P(^ e «-> z/ r ) 
o i/ T ) 



sin 2 c/>sin 2 2(9 sin 2 (A 3 iL/2) , 
cos 2 </>sin 2 26»sin 2 (A 3 iL/2) , 

sin 2 6 sin 2 20sin 2 (A 2 iL/2) + cos 2 6 sin 2 20sin 2 (A 32 L/2) 
-^sin 2 2(f) sin 2 26> sin 2 (A 3 iL/2) , 



(5) 
(6) 

(7) 
(8) 

(9) 



where A 3i = A 32 + A 2i . 

For large E, appropriate for atmospheric and long-baseline neutrinos, if a 2 <C (cE) 2 , then 
sin 2 9 <C 1 , cos 2 ~ 1 and the only appreciable oscillation is 



P(v^ o v T ) ~ sin 2 20sin 2 (A 32 L/2) 



where 



A 



32 



2cE 



(10) 



(11) 



Thus the — > u T oscillation amplitude has amplitude sin 2cj) and is maximal for <j) = f j m which 
case ai 2 = ai 3 (reproducing the simple two-parameter bicycle model). The energy dependence of 
the oscillation argument in this limit is the same as for conventional neutrino oscillations due to 
neutrino mass differences, with an effective mass-squared difference 

,2 



5m 



eff 



2EA 



32 



a 
c 



(12) 



The measured value for 5m 2 ^ in atmospheric and long-baseline experiments then places a constraint 
that relates a and c. 

If E is not too large, then the more general Eqs. (J3J)- dSD apply. Furthermore, in matter there 
is an additional term due to coherent forward scattering [13], which adds a \^2GpN e term to 
the upper left element of h e ff, where N e is the electron number density. In matter the angle (ft 
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is unchanged and 6 is now given by Eq. ([5]) with the substitution an — > an + GpN e /y/2. For 
adiabatic propagation in the sun the solar neutrino oscillation probability is 

P(v e u e ) = cos 2 6 cos 2 6»o + sin 2 9 sin 2 9 , (13) 

where 9q is the mixing angle at the creation point in the sun (with electron number density iVg ~ 
90 iV^/cm 3 ) and 6 is the mixing angle in vacuum. For convenience we define the quantity 6 = 
G F N°/{2V2) = 1.7 x 1(T 12 eV. 

The probability has a minimum value 

Pmin= 2^+62' (14) 

which is always less than ^. The minimum P must match the oscillation probability of the 8 B 
neutrinos (which from the SNO experiment |14j is P mm — 0.30), which fixes a to be 




2Pmin ) 

At very low energies the solar neutrino oscillation probability is 

1 



2Pm ' ! ' 2.1xlO- 12 eV. (15) 



low — 2 



an(axx + 2b) 



(16) 



^Ja 2 n + a 2 \/(aii + 2b) 2 + a 2 _ 

Note that the probability in Eq. (|16p is exactly | for an = (e.g., in the simple two-parameter 
bicycle model), which is not a good fit to the low-energy solar neutrino data. However, for an > 
or an < —26, the low-energy probability can be made larger than i. Using the low-energy value 
P ~ 1— | sin 2 2#i2 ~ 0.57, where 9\2 is the usual solar neutrino mixing angle [TS], we find an = 0.206 
or an = —2.26. 

The probability reaches the minimum at 

Emin = - [an + 6] , (17) 

which must occur in the energy region of the 8 B solar neutrinos (E m i n ~ 10 MeV), which fixes the 
magnitude of c to be 

1 ,,1-2 6 lq 
\c\ = - |on + 6| = - «2.0xl0- 19 . (18) 

Using Eq. (|12p we may now calculate the value of the atmospheric om 2 ^ inferred from solar neutrino 
data: 5m 2 = a 2 /c = 2.2 x 10 -5 eV 2 , which is two orders of magnitude below the measured value. 

One caveat for this calculation is that the low-energy solar oscillation probability is not measured 
precisely, and the model prediction may be adjusted by changing an. This in turn changes c (via 
Eq. [T8|) and the predicted atmospheric om 2 ^ (via Eq. [T2|) . The relationship between om 2 ^ and 
Plow is shown in Fig. [TJ where we have assumed 8 MeV < E m i n < 12 MeV and 0.27 < P m i n < 0.33. 
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Figure 1: Correlation of 5m 2 e ^ for high-energy atmospheric and long-baseline neutrinos with the 
oscillation probability Pi ow for low-energy solar neutrinos in the generalized direction-independent 
bicycle model (between the solid curves), found by varying the model parameter an. The left 
(right) solid curve assumes E min = 8.0 MeV (12 MeV) and P m in = 0.27 (0.33). The region between 
the horizontal dotted lines is consistent with atmospheric and long-baseline neutrino experiments. 

For the range of Sm^f allowed by experiment (shown by the horizontal dashed lines), the low- 
energy probability is approximately 0.30, which is not consistent with Pi ow « 0.57 preferred by the 
solar data. In fact, any Sm 2 ^ above 10~ 4 eV 2 gives a value for Pi ow below 0.40. Therefore there 
is no acceptable value for an that fits both the low-energy solar oscillation probability and Sm 2 ^ 
for high-energy atmospheric and long-baseline neutrinos, and the generalized direction-independent 
bicycle model is excluded. 

3 Other textures for h e ff 
3.1 Classification of models 

There are six possible c coefficients in h e ff. three real diagonal coefficients and three complex 
off-diagonal coefficients (the remaining three off-diagonals are fixed by the hermiticity of h e ff). 
Therefore there are 2 6 = 64 possible c textures for h e ff. Since the high-energy behavior of h e ff 
is determined by the c coefficients, we classify the models by the number of nonzero c there are 
in h e ff. Within each main class there are distinct subclasses which depend on the diagonal/off- 
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Table 1: A list of the 64 possible c textures for h e ff. The number in the subclass name corresponds 
to the number of nonzero c, while the letters indicate a distinct diagonal/off-diagonal structure (up 
to flavor permutation), if applicable. A Di in the structure column indicates that a diagonal ca is 
nonzero, while an Ojk indicates that off-diagonal Cjk is nonzero. Different latin indices in each case 
are distinct, e.g., in the structure DiOjk the diagonal element does not share a row or column with 
the off-diagonal element, whereas for DiOij it does. 



Number of 


Subclass 


Structure 


Number of flavor 


nonzero c 






permutations 








— 


1 


1 


1A 


A 


3 




IB 


On 


3 


2 


2A 




3 




2B 


DiOij 


6 




2C 


D l O jk 


3 




2D 


Oij Oik 


3 


3 


3A 


DiDjDk 


1 




3B 


1 );!),();, 


3 




3C 


DiDjOik 


6 




3D 


DiOij Oik 


3 




3E 


DjOijOik 


6 




3F 


OijOikOjk 


1 


4 


4A 


DiDjDkO^ 


3 




4B 


DiDjOijOik 


6 




4C 


DiDjOikOjk 


3 




4D 


DiOijOikOjk 


3 


5 


5A 


DiDjDkOijOik 


3 




5B 


DiDjOijOikOjk 


3 


6 


6 


DiDj DkOijOikOjk 


1 



diagonal structure; within each subclass there are textures that differ only by permutation of the 
flavor indices. In all there are 19 subclasses, which are listed in Table 1. 

We note that we may subtract any quantity proportional to the identity from h e ff, since common 
phases in the neutrino equations of motion do not affect the oscillations. In this way a diagonal 
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element may be removed or moved from one position to another. Then it is not hard to see that 
the following subclasses are strictly equivalent: 3Ao2A, 3C<H>3B, 4A<H>3B, 4Co4B, 5Ao4B and 
6o5B. 



3.2 Method for analyzing textures 

Our analysis proceeds as follows. We assume that \(%jE\ 3> \akt\ for any (i,j,k,£) for the high 
energies of atmospheric and long-baseline neutrinos. This assumption is justified since if any a 
is similar in magnitude to the cE at high energies, then at lower energies (such as for reactor 
neutrinos) the a terms will dominate and the oscillation arguments will be energy- independent, 
contrary to the KamLAND data, which measured a spectral distortion (similarly, solar neutrinos 
would also not have an energy-dependent oscillation probability, as they must). Furthermore, for 
the sake of naturalness, we assume that the c coefficients are all the same order of magnitude, and 
that likewise the a coefficients are also the same order of magnitude. 

Although for each texture the number of nonzero c is determined, initially we place no restric- 
tions on the a. We note that if all off-diagonal c are nonzero, then by a redefinition of neutrino 
phases and adding a term proportional to the identity we may take all off-diagonal c to be real and 
positive, except for one off-diagonal c that is complex (which we take to be C13 unless otherwise 
noted). If any off-diagonal c is zero, the nonzero off-diagonal c may all be taken as real and positive. 

A key feature of the bicycle model was that even though the terms in the effective Hamiltonian 
were either proportional to energy or constant in energy, one eigenvalue difference was proportional 
to E -1 , which mimics the energy dependence of the oscillations of atmospheric and long-baseline 
neutrinos. Having an eigenvalue difference proportional to E -1 means that if the eigenvalues are 
expanded in a power series in neutrino energy, 



then two eigenvalues must be degenerate at leading order in E (linear in E), and at the next 
order in energy (E°, independent of energy). Therefore in our analysis of more general three- 
neutrino models with Lorentz invariance violation, we look for model parameters that satisfy these 
conditions. Since an L/E dependence has been seen over many orders of magnitude in neutrino 
energy [16] . it seems likely that this is the only way the Hamiltonian in Eq. (prj will be able to fit 
all atmospheric and long-baseline neutrino data. 

For each texture we expand the eigenvalues of h e ff in powers of E (as in Eq. [T9|) . where the 
leading E 1 behavior comes from the dominant cE terms. Since we want 1/E behavior for at least 
one oscillation argument, we require that two of the eigenvalues be degenerate to order E°, with 
the first nonzero difference occurring at order E~ l . In all cases this requirement puts constraints on 



00 




for i = 1,2,3 



(19) 



3=0 



8 



the c and a coefficients. In our calculations we first find the eigenvalues to order E 1 and impose the 
constraint that two eigenvalues must be degenerate; then we find the eigenvalues of the simplified 
h e ff to order E° and again impose the degeneracy condition. In this way the expressions for the 
eigenvalues to order E~ x will be made as simple as possible at each stage of the calculation. 

If the appropriate 1/E behavior can be achieved, the mixing angles are then calculated to 
determine if f^s have maximal mixing and v e small mixing for atmospheric and long-baseline 
neutrinos. If the model is still viable, the energy dependences of the oscillations of solar and 
KamLAND neutrinos are then checked for consistency. 

At any time we are allowed to subtract a constant times the identity matrix from h e ff. Some 
cases may then be further simplified, or made equivalent to other cases (see below for specific 
examples). Rotations are also sometimes used to show that some cases are equivalent to others. 

3.3 No c parameters 

In this case, Class 0, h e ff has only a terms and therefore is independent of energy. This clearly 
cannot produce 1/E behavior at high energy, so this category is immediately ruled out. 

3.4 One c parameter 
3.4.1 Class 1A 

This case has the structure 



/ cE + an an «i3 ' 



h, 



eff 



12 



022 2 3 , (20) 



V °13 fl 23 °33. 

where c\\ = c may be taken as real and positive. The eigenvalues to order E° are then 

\i = cE + au, A 2 ,A 3 = i a 2 2 + a 33 ± \J (a 22 - G133) 2 + 4|a 2 3| 2 • (21) 

The difference A 2 — A3 can only be made zero to order E° if a 22 = 033 and 1 0,23 1 = 0. Then 033 
times the identity may be subtracted from h e ff, if an — 033 is redefined as an, this case reduces 
to the generalized bicycle model described in Sec. 2, which is excluded by the combined data. 

3.4.2 Class IB 

This case has the structure 

an cE + a\2 013 
h ef f = I cE + a* 2 a 22 a 23 | , (22) 

«*3 °23 a 33 . 
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where c\2 = c may be taken as real and positive. The eigenvalues to order E l are then 

Ai,A 2 = ±|c|£, A 3 = 0. (23) 

Since these are all different at leading order, they cannot give an oscillation argument proportional 
to E~ x at high energies, and this case is not allowed. 

3.5 Two c parameters 
3.5.1 Class 2A 

This case has the structure 



/ c n E + a u 



h ef.f = 



V 



12 



13 



ai2 an\ 

C22E + 0,22 O23 
033/ 



(24) 



A 23 



where c\\ and C22 are real. The eigenvalues at leading order are 

Xi = c u E, \2 = C2 2 E, A 3 = 0, 



(25) 



so that we must have cn = C22 for degeneracy (having one of the cu = also works, but then it is 
in Class 1A instead of 2A). Now if c\\E times the identity is subtracted from h e ff, this reduces to 
Class 1A, which is ruled out. 



3.5.2 Class 2B 

This case has the structure 

h, 



/cuE + an a 2 E + a 12 a n \ 



eff 



c 12 E + a\ 2 



\ 



13 



022 



A 23 



023 
(133/ 



(26) 



where c\\ and ci 2 may be taken as real and positive. The eigenvalues at leading order are 

1 



Ai, A2 — 



c u ± A/cfi +4c 



12 



E. 



A3 = . 



(27) 



Degeneracy requires (i) Ai = A2, which is not possible for nonzero c\\ and ci 2 , or (ii) A3 = Ai or 
A 2 , which is not possible for nonzero ci 2 . Therefore this case is not allowed. 

3.5.3 Class 2C 

This case has the structure 



h 



eff 



12 



>Cl3^ + a* 3 u, 23 

10 



ai2 c 13 E + ai 3 ' 

C22-E- 023 

a%i 033 



(28) 



where C22 and ci 3 may be taken as real and positive, and we have subtracted a term proportional 
to the identity so that 022 = 0. The eigenvalues at leading order are 



Ai, A2 = ^fcisE , A3 = C22E . 



(29) 



Degeneracy requires C13 = 022- If we define C22 = C13 = c, where c is a positive real number; then 
the eigenvalues to order E° are 



Ai = -cE + ^ [an + a 33 - 2Re(ai 3 )] 



A 2 ,A 3 = cE + 



1 



± yjx 2 + 8|y| 2 



(30) 



where x = an + 033 + 2Re(ai3) and y = au + a^- Degeneracy to order E° requires x = and 
y = 0, which implies an +033 = — 2Re(ai3) and a\2 = —023- With these conditions the eigenvalues 
to order E~ x are 

1 



Ai = -cE + an + 033 - 



2cE 



^2|a 2 3 1 + |ai3| -011^33), 



A, 



1 



cE + — ^-(2|a 23 | 2 + |ai 3 | 2 - a u a 33 ) , 
2cE 



A3 = cE . 



(31) 
(32) 



Clearly A32 = A3 — A2 has the correct energy dependence for atmospheric and long-baseline oscil- 
lations. The mixing matrix such that U T h e ffU is diagonal at leading order is given by 



U = 



1 


±sin0 


7S cos 







cos 


— sin( 


9 


1 

V2 


±sin0 


7S cos 





(33) 



where sin = \a\\ + 0-13 1 / \/ 2|a23 1 2 + |«n + «13| 2 and the oscillation probabilities are approximately 
given by 

1 



1 - sin 2 20 sin 2 ( -A 32 L 



i sin 2 20 sin 2 Qa 32 £ 



(34) 
(35) 



Therefore maximal mixing for is possible with Sm 2 ^ = 2E/S.23 = (2| 023 1 2 + |ai3| 2 — ana 33 )/c, 
but Up oscillates equally to v e and u T , which is excluded by atmospheric neutrino experiments. 
Hence this case is not allowed. 



3.5.4 Class 2D 

This case has the structure 



h, 



eff 



( an cnE + au c i3 E + a i3 \ 
C12-E + a\ 2 a 2 2 a 2 3 



\ci 3 E + 



(36) 



'13 



4 23 



fl33 
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where C12 and C13 may be taken as real and positive. If a rotation is applied to the fx — r sector, 
then ci 3 may be rotated away into ci 2 , which reduces this case to Class IB, which is not allowed. 
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3.6 Three c parameters 
3.6.1 Class 3A 

This subclass has nonzero c in each diagonal term and no off-diagonal c. By subtracting off c 33 E 
times the identity, this case reduces to Class 2A, which is not allowed. 



3.6.2 Class 3B 

This case has the structure 

h, 



/cuE + au ai2 ci 3 E + a n \ 



, I f - a ., a 23 , (37) 

\ c 13 E + a* 13 a% 3 c 33 E + a 33 / 

where en, C33 and C13 may be taken as real and 022 has been set to zero by a subtraction proportional 
to the identity. The eigenvalues at leading order are 

1 



Ai, A2 — 



Cu + C33 ± yj (c U - C 33 ) 2 + 4cf 3 



E. 



A3 = . 



(38) 



There are two possible ways to have a degeneracy. First, if Ai = A2, then we must have c\\ = C33 
and C13 = 0. However, if c\\E times the identity is then subtracted from h e ff, this possibility 
reduces to Class 1A. Second, we can have A2 = 0, so that it is degenerate with A3. There is a family 
of such solutions with C33 = r 2 c\\ and C13 = rcn, where r may be taken as a positive real number. 
If we define c\\ = c, then to order E° the eigenvalues are 

1 

2 



Ai = (1 + r 2 )cE + an + a 33 - x , 



A2, A3 



± \/x 2 + Ay 



(39) 



where x = [033 + r 2 a\\ — 2rRe(ai3)]/(l + r 2 ) and y = \rau — d2 3 \ 2 /(l + r 2 ). Degeneracy is only 
possible if x = and y = 0, which requires 033 + r 2 a\\ = 2rRe(ai3) and 023 = ra\ 2 ^ respectively. 
The eigenvalues to order E~ l are then 

|ai 3 | 2 + (1 + r 2 )\a 12 \ 2 - aua 33 



Ai = (1 + r z )cE + a u + a 33 + 



(1 + r 2 )cE 



A 2 = - 



N3I 2 + (1 + y 2 )|ai 2 | 2 - 011033 
(1 + r 2 )cE 

Thus A32 has the correct energy dependence, and gives 



A, 



0. 



Sm 2 eff 



2EA 



32 



2 |ai 3 | 2 + (1 + r 2 )|ai 2 | 2 - 011033 _ Jran - ai 3 | 2 + (1 + r 2 )|ai 2 | 2 



(40) 
(41) 

(42) 



(l + r 2 )c c(l + r 2 ) 

for atmospheric and long-baseline neutrinos. We note that A3 = is an exact result given the 
degeneracy conditions, true even when E is not large. 
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To leading order the mixing matrix that diagonalizes h e ff via U T h e ffU is 



U 



/ COS ( 





sin< 



8111 <p COS f 

sin# 

- cos <b cos i 




(43) 



where sine/; = r/yl + r 2 and tan# = \/T+ r2 |oi2|/l ra n — a i3|- This mixing gives the oscillation 
probabilities 



P(v L 



1 - sin 2 261 sin 2 | A 32 - 



P(v„ 



P(Ve 



sin 2 6 sin 2 26* sin 2 A S2 - 



(44) 
(45) 



L 



1 - cos^ 9 sin z 2^ shr A 2 i - - sin^ 6* shr 20 sin z A 3i - 



L 



sin 4 sin 2 2(9 sin 2 ( A 



L 



^32: 



(46) 



Maximal oscillations are possible for 9 ~ vr/4, which imposes the condition \/T+ ?~ 2 |ai2| — 
|ran - ai 3 |. 

Oscillations of z/ e at high energies must be small due to the limit on v„ —> v e from K2K [T7] 
and MINOS [18J0 For K2K and MINOS the oscillation amplitude for P(u^ -> u e ), sin 2 0sin 2 29, 
has an upper bound of about 0.14, which implies r < 0.43 for 9 ~ ir/4. The T2K experiment sees 
evidence for Un —> v e at the 2.5a level |19| : their allowed regions are consistent with this bound. 

We note that the conditions C33 = r 2 cn and C13 = rc\\ require fine tuning. If these conditions 
are not exact, they introduce small corrections, which may be absorbed into the a terms, e.g., 
dij — > a,ij + 5cijE, where 5cij represents the deviation from the exact degeneracy condition. This 
effectively introduces an E dependence into brr? e ^, contrary to the atmospheric and long-baseline 
data. 

For solar or reactor neutrinos the large energy limit does not apply. Then the eigenvalues are 
l r 



Ai, A 2 
A3 



2 

0. 



cE(l + r 2 ) + a u + a 33 ± J[cE(l + r 2 ) + a u + a 33 ] 2 + 25m 2 //C (l + r 2 



(47) 



where Sm^rr is from Eq. (|42p . and it can be shown that the matrix that diagonalizes h e ff is 



/ cos cos £ + sin cos 9 sin £ e~ tS 
U = sin 9 sin £ 

\ sin cos £ — cos cos 9 sin £ e _i<5 



cos sin £ + sin cos cos £e~ J<5 — sin sin 9 \ 

sin cos £ cos 6>e* 5 

sin sin £ — cos cos # cos £e~~ J<5 cos sin / 



, (48) 



1 Limits on — > !/ e or ^ M — >• £> e from experiments such as CHOOZ or KARMEN do not apply here since they 
involve lower energy neutrinos. MiniBooNE limits may apply, but only for Sm^ff <; 0.1 eV 2 , and therefore not at 
the A32 scale. 
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Figure 2: The best fit to the KamLAND data for Class 3B (solid lines) and the standard oscillation 
scenario with neutrino masses (dashed lines). For the model parameters see Eq. (|49p . 

where (p and 9 are defined as above, tan£ = \ra\% — ai3|/(Ai cos#) and 5 = arg(ran — an). Note 
that in the large energy limit Ai is large, £ — > 0, and Eq. (|48p reduces to Eq. (|43p . Also, since none 
of the mixings are zero, CP violation is possible. 

We checked KamLAND phenomenology first. Since we have obtained several conditions from 
fitting the atmospheric and long-baseline neutrinos, using these conditions we can vary an, ai3 and 
r to fit the KamLAND data [20]. Other parameters in the effective Hamiltonian will be determined 
by these three parameters. Scanning the an, ai3 and r parameter space, we find the following 
parameter values yield reasonable agreement with the KamLAND data (see Fig. [2]): 

a n = -8.7 x KT 11 eV , o 13 = -9.5 x KT 11 eV , r = 0.1 . (49) 

However, the fit is not as good as the standard oscillation scenario with neutrino mass. 

Next we use these parameter values to check the solar phenomenology. Since the operator for a 
breaks CPT, we reverse the sign of a when we apply these parameter values to the solar neutrinos. 
However, the prediction does not agree with the solar data at high energies given the upper bound 
on r from above (see Fig. [3|). 

We also searched the an, a\3 and r parameter space to fit the solar data separately. The best 
fit still can not yield reasonable agreement with the solar data at high energies for r < 0.43 (see 
Fig. If we do not impose the constraint on r, the fit to the solar data is improved at high 

In order to understand why the oscillation probability for the high-energy solar neutrinos is so high, we consider 
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Figure 3: The prediction of Class 3B for the oscillation probability of solar neutrinos using the 
parameter values obtained from fitting KamLAND data [20J. The solar data points are from an 
update of the analysis in Ref. [21]. 



energies (see Fig. [5]). However, we cannot simultaneously fit the KamLAND and solar data even 
with larger r. We found that we also need |on| to become larger in order to fit the solar data, but 
larger \a\i\ yields fast oscillations for KamLAND data with averaged probabilities around 1/2. 



the survival probability of solar neutrinos in the high energy limit. As we have noted, the mixing matrix in vacuum 
reduces to Eq. (I43[) in the high energy limit. In matter, we can still write the mixing matrix in the form, 

/ cos <j) sin cos #o — sin <f> sin 8q ' 



U = 



sin 6*o cos 6q 

\ sin <f> — cos 4> cos do cos <f> sin Go 



(50) 



except tan So = ry/l + r 2 |ai2|/|a33 — ra*3[, since we do not have the relation 033 + r 2 a\\ — 2rRe(ai3) in matter (but 
0q is still the same as 6 in vacuum). Now the survival probability of the solar neutrinos in the high energy limit is 



Since 80 = 9 ~ n/4 in vacuum, we have 



cos 4 4> + — sin 4 (f>{l + cos 28 cos 28q) 



4 1 1 1 -4 

cos sm ' 



3, . 2 
-(sm , 



3' 



(51) 



(52) 



Since sin</> = r/y/1 + r 2 , applying the constraint for r, r < 0.43 gives sin 2 < 0.14, and the v e survival probability 
approaches 0.75 in the high energy limit. This is the reason that we cannot fit the solar data at high energies with 
the constraint on r. 
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Figure 4: Best fit prediction for survival probability of v e for solar neutrinos for Class 3B, assuming 
r < 0.43. The model parameters for the best fit are an = —1.0 x 10~ 10 eV, a i3 = 7.3 x 10 -11 eV 
and r = 0.4. 

1.0 i 1 




0.0 I ' ■ — ■ ■ — 1 

0.1 1 10 

E v {MeV) 

Figure 5: Best fit prediction for survival probability of v e for solar neutrinos alone in Class 3B, 
assuming r > 0.43. The model parameters are an = —5.0 x 10~ 10 eV,ai3 = —3.1 x 10 -11 eV and 
r = 1.0. 
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3.6.3 Class 3C 



This case has the structure 



/c n E + an cti2 c 13 E + a l3 ' 



h, 



eff 



b 12 



C22E + 022 0-23 

\a 3 E + a* 13 a* 23 a 33 

where c\\ and C22 are real and c± 3 may be taken as real and positive. By subtracting C22 times the 
identity, this case reduces to Class 3B, which was described in the previous section. 

3.6.4 Class 3D 

This case has the structure 



(53) 



h, 



eff 



fcuE + au cnE + au ci 3 E + a 13 \ 
c 12 E + a\ 2 a 2 2 a 2 3 



\ci 3 E + 



(54) 



13 



4 23 



033 



where en, c\2 and c\ 3 may be taken as real and positive. If a rotation is applied to the ji — r sector, 
then C13 may be rotated away, which reduces this case to Class 2B, which is not allowed. 

3.6.5 Class 3E 

This case has the structure 

/ cnE + a u ci 2 E + a 12 ai 3 \ 
h e ff = cuE + a* 12 a 2 2 c 23 E + a 2 3 , (55) 

\ a* 13 C2 3 E + a* 23 a 33 / 

where en, C12 and C23 may be taken as real and positive. This is the first case that cannot be simply 
reduced to a previous case, and which requires solving a nontrivial cubic equation to determine the 
eigenvalues at leading order. The eigenvalue equation for h e ff/E at leading order is 



A 3 - ciiA 2 - (c\ 2 + c\ 3 )\ + c u c 2 23 = . 
For a cubic equation of the form A 3 + a A 2 + 6A + c = 0, if we define q 



(56) 



36 and r 



2a 3 — 9ab + 27c, then for three real roots the cubic discriminant / = 4g 3 — r 2 must be nonnegative, 
with / = when two of the roots are equal. Since the effective Hamiltonian is hermitian, the 
eigenvalues must be real, so / > 0. Therefore, if there is a degeneracy, not only must / = 0, it 
must be a global minimum of /, i.e., we can search for degeneracies by finding the minima of /. 
For this case we have 



9 = 4 + 3(c 2 2 + 4 3 ) , 



r = cn(-2c? 1 + 18^3 



9c? 2 ) , 



(57) 
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and the discriminant is 



/ = 4(4 + 3c 12 + 34) 3 - c 2 u (2c?! + 9cj 2 - 18c 2 23 ) 2 . (58) 

Then, 

= ^- = 108c 12 [64 + 64 + 44 + 10c?! + 1244] > (59) 
requires at least that C12 = 0, which reduces this case to Class 2C, which is not allowed. 

3.6.6 Class 3F 

This case has the structure 



h eff 



/ «n c 12 E + a 12 ci 3 E + a 13 ' 

ci 2 E + a* 12 a 22 c 23 E + a 23 \ , (60) 
\cl 3 E + a* 13 c 23 E + a* 23 

where c\ 2 and c 23 may be taken as real and positive, c\ 3 is complex and a 33 has been set equal to 
zero. This case also requires solving a nontrivial cubic equation to find the eigenvalues at leading 
order; with 

Q = 3 ( c i2 + c i3 + c 23) > r = -54ci 2 c 23 |ci 3 |c 5 , (61) 
where c$ = cos 5 and 8 is the phase of c\ 3 . Searching for a minimum of / = 4g 3 — r 2 : 

= -^- = 72ci2<? 2 + 108r|ci3|c23C5, (62) 



Oc 



12 



= ^ = 72|ci3k 2 + 108rci2C 23 c 5 , (63) 
dc 13 

df 

= ^ = 72c 23 q 2 + W8rc 12 \c 13 \c s , (64) 

OC23 
df 

= — = -108rci2|ci 3 |c23sm5. (65) 

The quantity q is explicitly nonzero; if r was zero then Eqs. (|62p - (|64p would imply that c\ 2 , c\ 3 and 
c 23 would all have to be zero, which is not possible for this case, so r ^ 0. Then the last equation 
implies sin<5 = 0, or 5 = or ir. Thus c$ = ±1, i.e., c\ 3 is real, but it might differ by a sign from 
C12 and C23; we use c$ to denote this possible sign difference and henceforth take C13 as real and 
positive. 

It is not hard to show that C12 = c\ 3 = c 23 is required for a minimum of / and that this condition 
gives / = 0. Therefore degeneracy requires c\ 2 = c\ 3 = c 23 = c with cs = ±1. Then the eigenvalues 
to order E° are 

l r 



Ai = 2cEcs + an + a 22 - x , A 2 , A3 = -cEcs + - 



=F Vx 2 - 4y , (66) 
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where 

2 

x = - [(a u + 022) - Re(ai 3 ) - c 5 Re(ai 2 + a 23 )] , , (67) 
V = g [ana22 - |ai2| 2 - |ai3| 2 - |a 2 3| 2 + 2Re(ai2«23 - 022013) 

+ 2c 5 Re(ai 2 ai3 + a 2 3a* 13 ~ «n«23)] • (68) 

Thus degeneracy requires that the quadratic discriminant g = x 2 — 4y be zero. Since the eigenvalues 
are real, we know g > 0, and degeneracy can only occur at a minimum of g. It can be shown that 
g has a minimum at zero for an = 022 = Re(ai2) = Re(ai3) = Re(a23) = and Im(ai3) = 
C5lm(ai2 + 023)- Then 

/ cE + ia\2 cs[cE + i(a i2 + 023)] \ 

h e ff = cE — iau cE + 2023 , (69) 

\c s [cE - z(ai 2 + 023)] cE-ia 2 3 / 

where the a^- are now defined as real. The eigenvalues of this matrix to order E' 1 are 

2cs 2c$ 
Ai = 2cscE+-^-^(aj 2 + al 3 +a 12 a23) , A 2 = -c s cE , A 3 = -c s cE- -^-^(a 2 12 + al 3 + ai 2 a23) , 

(70) 

and the mixing matrix that diagonalizes h e ft is 

U = V§ -%( a u + a ™) 7^(012-023) , (71) 

where N = \f2(a\ 2 + O23 + 012O23) is a normalization factor. 

At high energies for the fast oscillation A31 ~ A21 — —3c$cE, all off-diagonal oscillation prob- 
abilities have the same approximate form: 

n 4 2 (3cEL\ 
P( 1/a ^^) = _ sin 2/__j . (72 ) 
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For this oscillation amplitude, 4/9, the NuTeV limit on — > v e [22J gives a 90% C.L. upper bound 
on Sm 2 ff of 3.6 eV 2 . Since we have ^m^ff = 6cE 2 in this case and the average neutrino energy was 
74 GeV, the experiment imposes the upper bound c < 1.1 x 10 -22 □ On the other hand, in order 
for the expansion in powers of E to be valid, we need N/ (cE) <C 1 for E > 100 MeV, which leads 
to the lower bound c > 4 x 10~ 17 . Therefore the structure required for the 1/E behavior at high 
energy is inconsistent with accelerator bounds. Since all flavors have the same survival probability 
in the fast oscillation, the result is the same even if a different permutation of flavors is taken. 



3 The NuTeV bound on Sm 2 is not the most stringent for ordinary massive neutrino oscillations, but because 
8m 2 e ff oc E 2 , the high neutrino energies in NuTeV give the strongest bound on c. 
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Furthermore, for the A23 oscillations in atmospheric and long-baseline neutrinos, all three flavors 
have the probability 

P(»a = ^ a2 \ 2 (| " l^l 2 ) ^ (^p) • (73) 

This implies that all flavors of downward atmospheric neutrinos would be suppressed by a factor 
of 5/9, which is contrary to the data. Therefore this case is excluded. 

3.7 Four c parameters 

3.7.1 Class 4A 

This case has three nonzero diagonal and one nonzero diagonal c. By subtracting a piece propor- 
tional to the identity, this case may be reduced to either Class 3B or 3C. 

3.7.2 Class 4B 

This case has the structure 



/cuE + an c\ 2 E + a 12 c 13 E + a 13 \ 



h 



eff 



(74) 



C12-E + a* l2 C22E + a 2 2 a 2 3 
\c 13 E + a* 13 a* 3 a 33 / 

where en and C22 are real and ci 2 and C13 may be taken as positive. The eigenvalue equation for 
h e ff/E at leading order is 

A 3 - (en + c 22 )A 2 + (ciic 22 - c? 2 - C2 3 )A + c 22 c 2 3 = . (75) 

This case has cubic discriminant 

f = 4q 3 -r 2 , (76) 

where 

<? = Cn + c| 2 - C11C22 + 3ci2 + 3c 2 3 , (77) 

r = -2cf 1 + 3c? 1 C22 + 3ciicl2-2c^2-9cii(c?2 + c? 3 ) + 9c22(2c?3-c?2), (78) 

The minimum conditions are 

= ^- = 12(2cn - c 22 )g 2 - 2r [-6c 2 ! + 6ciic 2 2 + 3c 2 2 - 9(c 2 2 + c 2 3 )] , (79) 

= = 12(2c 22 - c n )q 2 - 2r [-6c 2 2 + 6c 22 cii + 3c 2 u + 9(2c 2 3 - c 2 2 )] , (80) 

CC22 

= £L = 72c 12 q 2 + 36rc 12 (c n + c 22 ) , (81) 



dc 12 
dc 13 



= ^ = 72c 13 g 2 - 36rci3(2c2 2 -cii). (82) 
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Clearly q > if none of the are zero. If r = 0, then Eqs. (|81j) and ([82]) would imply C12 = C13 = 0, 
which is not Class 4B; therefore r ^ 0. Then Eqs. (18ip and (|82p imply 

1 r 1 



en + c 22 q 2 2c 22 - en ' 
which implies C22 = 0. This case then reduces to Class 3D, which is not allowed. 

3.7.3 Class 4C 

This case has the structure 



(83) 



h eff 



/c n E + a n a 12 c 13 E + a 13 

a* 2 C 22 E + a 22 023-^ + 023 | , (84) 

\c 13 E + a* 13 c 23 E + a* 23 a 33 
where c\\ and C22 are real and c\ 3 and C23 may be taken as real and positive. By subtracting en 
times the identity this case may be reduced to Class 4B, which is not allowed. 

3.7.4 Class 4D 

This case has the structure 

I cnE + an c\ 2 E + a 12 ci 3 E + a 13 \ 
Kff = c X2 E + a\ 2 a 22 c 23 E + a 23 , (85) 

\c* 13 E + a* 13 c 23 E + a* 23 a 33 J 
where en, c\ 2 and C23 may be taken as real, and c\ 3 is complex. The eigenvalue equation for h e ff/E 
at leading order is 

A 3 - cnA 2 + ~{c\ 2 + c\ 3 + c 2 23 )\ + C11C23 - 2c 12 ci 3 c 23 c s = , (86) 

where cs = cos 5, c\ 3 — > c\ 3 e %& and c\ 3 is now taken as real and positive. This case has 

q = c 2 n +3(c 2 12 +c 2 13 + 4 3 ), (87) 
r = -2c 3 u + 9c u (2cl 3 - c 2 12 - c 2 13 ) - 54ci 2 c 13 c 23 c s , (88) 

where the discriminant is / = 4g 3 — r 2 . The minimum conditions are 

= ^L= 24c u q 2 ~ 2r [-6c 2 n + 9(2c 2 23 - c\ 2 - c 2 3 )] , (89) 

= - — = 72c 12 q 2 + 36r(cnci2 + 30130230,5) , (90) 
oci 2 

Qf 

= - — = 72ci 3 q 2 + 36r(dici 3 + 3ci2C2 3 c 5 ) , (91) 

CC13 
Q f 

= - — = 72c 23 q 2 - 36r(2ciic 2 3 - 3ci 2 ci 3 c s ) , (92) 
CC23 

df 

= — = -108rci 2 c 13 c 23 sin5. (93) 
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Clearly q > if none of the c are zero. If r = 0, then Eqs. (|90|) - (|92|) would imply c\i = C13 = C23 = 0, 
which is not Class 4D; therefore r / 0. Thus Eq. (193j) implies sin<5 = 0, or cs = ±1; therefore the 
off-diagonal elements are real. 

By combining Eqs. (f90|) and (f9Tj) . we find C12 = C13, and by combining Eqs. (f90|) and ([92]) . we 
find c 2 2 = c 2 3 + C11C23C5. Then g = (en + 3csC23) 2 and r = — 2(cn + 3C5C23) 3 . Clearly then / = 0, 
and the conditions for degeneracy at leading order are 



C12 = C13 



c 12 = c 23 + c <5CllC 2 3 



(94) 



Thus there is a two-parameter set of degeneracies at leading order for this texture; at leading order 
h e ff has the form 

/ en S Sc s ' 



h 



eff 



S c 23 I E , 

\Sc S C 2 3 



(95) 



where S = £23(023 + c^Cn). By applying the rotation 



V 



(I 



(96) 



and adding a term CSC23E times the identity, at leading order the new Hamiltonian is 

/C11 + CSC23 V2Scs' 



hi 



eff 



V T heffV 



I E. 

\ V2Sc 8 2c 5 c 2 3 



(97) 



Equation (|97p has the form of Class 3B with r = \j2cizl {c-iz + QCn). The matrix that diagonalizes 
the original h e ff is therefore U' = VU, or 

/V2cos(^ \/2 sin <f> cos 9 — 72sin</>sin# \ 

U' = —= cs sin <p sin 9 — cs cos cos 9 cos 6* + C5 cos 4> sin , (98) 
\ cos (ft — cs sin — cos (ft cos 9 — cs cos 9 + cos 4>sm9 ) 

where U is from Eq. (|43[) . The oscillation probabilities are 



1 — (sin# — cos </>c,5 cos 6>) z (cos 9 + cos </>ca sin#) z sin^ ( A23— - 



sin 2 (j> (sin 6 — cos cos 6>) 2 sin 2 ( A12 — - 



L 
2 
L 
2 

sin 2 sin 29 (sin 2 sin cos 9 — cs cos </> cos 29) sin 2 I A23- 



sin 2 cf> (cos 9 — cs cos </> sin 0) 2 sin 2 ( A13 



(99) 
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— eg sin 4> sin 2(/>(sin 9 cos 9 — eg cos <\> cos 2 6*) sin 2 ^Ai2 — 
+c,5 sin (/>sin 2(/>(sin6' cos + Q cos 0sin 2 0) sin 2 ^Ai3— ^ , (100) 
P(i/ e ->i/ e ) = 1 - sin 4 0sin 2 26»sin 2 ^A 23 ^ - sin 2 20cos 2 6»sin 2 ^Ai2^ 

- sin 2 20 sin 2 # sin 2 (&13^) • (!01) 

In order to compare with the atmospheric and long-baseline neutrinos data, for large E, we 
should have A12L, A13L 3> A 2 3.L ~ 1. Then the oscillation probabilities are 

P(Vfi, — > Vy) = 1 — (sin# — eg cos 0cos 9) 2 (cos 9 + eg cos 4>sm9) 2 sin 2 ^A23 — 

sin 2 0(1 + cos 2 <p), (102) 



P{Vfj, v e ) = sin 2 (f> sin 2$(sin 2 sin 9 cos — eg cos cos 20) sin 2 ^ A23— ^ +- sin 2 2o .( I ()■'!) 
Maximal f e oscillations requires 



1 ~ (sin 9 — eg cos cos 9) (cos + c<j cos sin 0) 

1 1 

= 1 - sin 2 4> - (cos 2 + -eg cos - - sin 4 0) sin 2 26> (104) 
If is small and sin 29 ~ 0, the probabilities are appropriate for the atmospheric and long- 



baseline neutrinos. Since sin0 = r f\J\ + r 2 and tan# = vT + r^a^l/lran — 013) > this imposes the 
conditions: (i) r ~ and a\2 — 0, or (ii) r ~ and ran — a i3- 

Since this case is equivalent to Class 3B after a rotation in the — v T sector, the v e — > v e 
oscillation probability expression is still the same. The results are also similar to Class 3B. While 
there are parameter values that yield reasonable agreement with the KamLAND data (see Fig. [6]) , 
they did not agree with the solar data at high energies (see Fig. [7]). 

Also, we fit the solar data separately. As was the case for Class 3B, we do not find a good fit 
to the solar data at high energies (see Fig. [8]) . 

3.8 Five c parameters 
3.8.1 Class 5A 

This case has three diagonal and two off-diagonal nonzero c. By subtracting a piece proportional 
to the identity, this case may be reduced to 4B or 4C. 
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Figure 6: Class 4D (solid lines) and the standard scenario (dashed lines) compared to the Kam- 
LAND data. The model parameters are on = —9.0 x 10 -11 eV, a±2 = 0, ars = —9 x 10 -11 eV and 
r = 0.02. 
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Figure 7: The best fit prediction of Class 4D to the solar data. The model parameters are an = 
9.0 x lO" 11 eV, ai2 = 0, a 13 = 9 x lO" 11 eV and r = 0.02. 
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Figure 8: Best fit prediction for survival probability of solar v e for Class 4D. The model parameters 
are a n = -6.0 x KT 11 eV, a 12 = 7.0 x 10~ n eV, a 13 = 0.42 x lO -11 eV and r 



-0.07. 



3.8.2 Class 5B 

This case has the structure 

h eff = 



fcuE + au cuE + au ci 3 E + a 13 \ 

c\ 2 E + a\ 2 a 22 c 23 E + a 23 
V c 13 E + a\ z c 23 E + 4 3 c 33 E + a 33 J 



(105) 



where c\\ and c 33 are real, c\ 3 and c 23 may be taken as real and positive, and c\ 2 is complex. At 
leading order the cubic equation for the eigenvalues of h e ff/E is 



A 3 - A 2 (cn + c 33 ) + A(cnc 33 - c\ 2 - c( 3 - C23) + C11C23 + c 22 cj 3 - 2ci 2 ci 3 c 23 c s , 

where the c\ 2 is the magnitude and 5 the phase of c\ 2 . Then we have 

q = c 2 n + c 2 33 - C11C33 + 3(c? 2 + c? 3 + c 2 23 ) , 
r = -2(cn + c 33 ) 3 + 9ciic 33 (cii + c 33 ) + 9c n (2c23 - cf 2 - ci 3 ) 
+9c 33 (2ci2 - cj 3 - C23) - 54ci 2 c 13 c 23 c 5 . 

and the conditions for a minimum of / = 4g 3 — r 2 are 

= ^ = 12(2c 11 -c 3 3)(7 2 -2r[-6c? 1 + 6c llC 33 + 34 + 9(2ci3- C ?2-c?3)] 
Of 

= 7 ^ = 12(2c 33 -c u )q 2 -2r[-6c 2 33 + 6c 11 c 33 + 3c 2 11 + 9(2c 2 12 -c 2 3 -c 2 23 )] 
CC33 



(106) 



(107) 
(108) 
(109) 



(110) 
(111) 
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df 

= ~pT~ = Vcm 2 + 36r [(en + c 33 )c 13 + 30120230,5] , (112) 

CC13 

= - — = 72ci2<? 2 + 36r [(en - 2c 33 )ci2 + 3ci 3 c 23 c,5] , (113) 

OC\2 
Qf 

= -z — = 72c 23 g 2 - 36r [(2cn - c 33 )c2 3 - 3ci 2 ci 3 c 5 ] , (114) 
wc 23 

df 

= — = -108rci2Ci 3 c 23 siii(5. (115) 
00 

It can be shown by the usual arguments that q and r are not zero, and cs = dbl. By eliminating 
g and r from Eqs. f 1 11 2 () and (|lX3j) we find 

C 33 Cl 2 Cl 3 = C 5 C 23 (Ci 3 - C? 2 ) , (116) 

and applying a simlar procedure to Eqs. (|1 12j) and ([1140 gives 

cnci 3 c 23 = c 5 ci 2 (ci 3 - C2 3 ) . (117) 

Using the relations in Eqs. (|116|) and (jll7j) it can be shown that all minimum conditions are met and 
that / = 0, so Eqs. ([1160 and ([1170 are the degeneracy conditions. Thus there is a three-parameter 
set of degeneracies at leading order for this texture. Then, after adding the term C12C23C5 / 'c\ 3 times 
the identity, the effective Hamiltonian at leading order may be written as 



h eff = 



Cl2Cl 3 C 2 3 



/ c 12 c 13 C13C23C12 C5C? 3 Ci2C23 \ 

Cl 3 C 23 cf 2 cf 2 C^ 3 C 5 Ci2Ci 3 C^ 3 E . (118) 

\cscl 3 C 12 C23 C6C 12 C 13 C% 3 cfgcfg / 

Without loss of generality 022 may be set equal to zero. Then the eigenvalues of h e f t to order E° 
are 

± v^ 2 - 4y] , (119) 



Ai — E + an + 033 1 ^2, A 3 — - 

Cl2Cl 3 C 23 2 



where 



- [anc^(c 2 2 + c 2 3 ) + a 33 c 2 2 (c 2 3 + c§g) 

-2Re(ai2)c 2 2 ci 3 c 23 - 2Re(ai 3 )c 2 3 ci2C 23 c 5 - 2Re(a 23 )ci 2 ci 3 C2 3 c 5 ] , (120) 

y = 02 [- 2a nRe(a23)ci2Ci3C2 3 C5 - 2a 33 Re(ai2)ci 3 C23C 2 2 + aiia 33 c 2 2 C2 3 

- c i3 c 23l a i2| 2 - c? 2 c? 3 1 a 23 1 2 - c 2 2 C2 3 |ai 3 | 2 + 2Re(ai 3 a 23 )ci 3 C2 3 c? 2 

+2Re(ai 3 ai 2 )ci 2 ci 3 C2 3 c 5 + 2Re(a2 3 ai 2 )ci2C2 3 Ci 3 c (S ] , (121) 

and S = c\ 2 c\ 3 + c\ 2 c% 3 + cf 3 c| 3 . Thus degeneracy requires that the quadratic discriminant g = 
x 2 — Ay be zero. It can be shown that g has a minimum at zero when 

( 2/13 2/23 A, / 100 n 

2/12 = c 5 ci2 , (122) 

V C13 C 2 3 / 
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C23011 . Ci 2 (cf 3 + C^gJ / C5C12C33 



112 = l^ + ^^ + cy (123) 

2C12C 23 C 23 (cf 2 +Cf 3 ) V 2c 13 / 

where = Re(ojj), and g/jj = Im(ajj); these are the degeneracy conditions. The eigenvalues to 
order E~ l are then 

c SE 1 

Ai = 1- -777^0,5012013023 [011033 - |«12 1 2 - |ai3| 2 - |a23| 2 l , (125) 

C12C13C23 SE 

A2 = 77g c <5Ci2Ci3C23 [011033 - |ai2| 2 - |ai3| 2 - |a-23| 2 ] , A 3 = , (126) 

and to leading order the mixing matrix that diagonalizes h e ff is 

/ cos 9 cos (ft — cos £ sin (ft — sin £ sin 9 cos sin £ sin — cos £ sin cos (ft \ 
{7= sin# cos # sin £ cos cos £ , (127) 

\ cos # sin (ft cos £ cos — sin £ sin sin (ft — sin £ cos (ft — cos £ sin 9 sin / 



where 



sin6> = -3=ci 2 C23, sin0 = —= ^ 23 _ , (128) 



V 5 V c 12 + c 23 



and 



sine S c /^3(an + a33) ) ^an-^ 

^3^2+4* ^ ^23 + C 2 2 



^3 = a 33 c 12( c 13 + c 23) + «ll c 23( c 12 + c 13) ~ 2aiia33C 13 C 12 C 2 3 , (130) 



is a normalization factor. The oscillation probabilities are 

P{v^^y^) = 1 - sin 2 £ sin 2 2(9 sin 2 ^A 2 i^ - cos 2 £ sin 2 26* sin 2 ^A 3i ^ 

- cos 4 9 sin 2 2£ sin 2 (a 32 ^ , (131) 

P(y^ — * v e) = 2 sin 29 cos 9 cos (ft sin £(cos £ sin (ft + sin £ sin cos </>) sin 2 ^A2i — 

—2 sin 29 cos 6* cos cos £ (sin £ sin (ft — cos £ sin 9 cos 0) sin 2 ^A3i — 
+2 sin 2£ cos 2 #(cos £ sin </> + sin £ sin cos (ft) (sin £ sin — cos £ sin cos (ft) 
x sin 2 (a 32 §J • (132) 

In order to have nearly maximal oscillations at the atmospheric scale, the A32 term must 
have amplitude close to unity, or 9 ~ 0,7r and £ ~ ir/4. Then from Eq. (|129p 

a 3 3Ci2(ci2 + C5C 2 3sin6') ~ aiic 2 3(c 2 3 - ^012 sin 9) . (133) 
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The small value of sin 9 implies cf 2 ,C2 3 <C c 2 3 . Furthermore, in order to have small Un — > v e 
oscillations at the A32 scale, sin 2 <p <C 1, or <C cf 2 , i-e., there is a hierarchy among the off- 
diagonal Cij. Then Eq. (|133p implies 033 <C an as well. Therefore there is a lot of fine tuning 
required to achieve the proper mixing. 

For simplicity, we only considered the parameters as real numbers. We have scanned the C12, C13 
and C23 parameter space to fit the KamLAND and solar data. Other parameters in the Hamiltonian 
can be determined by these three parameters, i.e., c\\ and C13 can be determined from Eqs. (I116P 
and (I117p . Also, for the atmospheric and long-baseline neutrinos, the A23 term has the correct 
energy dependence, and gives 

2 

Sm eff = 2EA 23 = ^C,5Ci2Ci3C23[aiia33 ~ 1^12 1 2 ~ |<3i 3 1 2 - | a 2 3 1 2 ] • (134) 

The above equation together with Eqs. (|124[) and (|133|) determine all a^. Another constraint is 
the hierarchy among the off-diagonal Cy, <C cf 2 *C cf 3 , which is also considered during the 
parameter search. 

We have varied the range of C13 from the order of 10™ 20 to 10~ 16 and take c\i and C23 to be at 
least one order of magnitude less than C13 and c\2 respectively. We found parameter values that 
can fit the KamLAND data (see Fig. [9]) , but they do not yield reasonable agreement with the solar 
data at high energies (see Fig. I10p . We also attempted to fit solar neutrinos alone and found there 
are no parameter values that can yield reasonable agreement with the solar data. The best fit is 
shown in Fig. [TT1 

3.9 Six c parameters 

In this case all c elements are nonzero. By subtracting off a piece proportional to the identity, this 
case may be reduced to Class 5B, which is ruled out. 

4 Summary 

We have examined the general three neutrino effective Hamiltonian in Eq. ([1]) for the case of 
direction-independent interactions and no neutrino mass. We looked for texture classes in which two 
eigenvalues were degenerate to order 1/E at high neutrino energy, so that oscillations of atmospheric 
and long-baseline neutrinos would exhibit the usual L/E dependence. 

Among the classes that had the proper 1/E dependence at high energy, none was also able 
to fit the atmospheric, long-baseline, solar and KamLAND data simultaneously. Class 1A (along 
with the equivalent Classes 2A and 3A) reduced to the direction-independent bicycle model, which 
has been shown to be inconsistent with the solar, atmospheric and long-baseline neutrino data. 
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Figure 9: Best fits for Class 5B (solid lines) and for standard oscillations (dashed lines) compared 
to the KamLAND data. The model parameters are C23 = 5.9 x 10~ 23 , C12 = 1.0 x 10 -22 and 
C13 = 8.9 x KT 19 . 
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Figure 10: The prediction of Class 5B for the solar neutrino survival probability using the parameter 
values obtained from fitting the KamLAND data. 
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Figure 11: Best fit prediction for the v e survival probability in Class 5B. The model parameters for 
the best fit are c 23 = 1.0 x 10~ 18 , c 12 = 2.4 x 10~ 18 and c 13 = 1.0 x 10" 18 . 

Classes 2C (and the equivalent 3E) and 3F did not have the proper oscillation amplitudes for 
atmospheric neutrinos. Finally, Classes 3B (and the equivalent Classes 3C, 4A and 4D) and 5B 
(and the equivalent Class 6) were able to fit atmospheric and long-baseline neutrino data, but could 
not simultaneously fit KamLAND and solar data at lower neutrino energies. The major difficulty 
in these latter classes was reproducing the low survival probability of high-energy solar neutrinos. 

Although we have not made an exhaustive search of the parameter space, the fact that high- 
energy neutrinos exhibit an L/E dependence in their oscillations over many orders of magnitude in 
E suggests that the only way this can occur in the effective Hamiltonian described by Eq. ([1]) is via 
the degeneracy of two eigenvalues to order 1/E. Since none of the cases where such a degeneracy 
occurs are also able to fit all neutrino data simultaneously, it seems extremely unlikely that any 
direction-independent SME model without neutrino mass will provide a viable description of all 
neutrino oscillation phenmomena. There is also strong evidence against direction-dependent terms. 
Furthermore, nonrenormalizable Lorentz noninvariant effective Hamiltonians with higher powers 
of energy (as in, e.g., the model of Ref. [9]) and no neutrino masses would require additional 
degeneracy conditions. Therefore it appears highly unlikely that Lorentz invariance violation alone 
can account for all of the observed oscillation phenomena. 
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